Complex-g' zeros of the partition function of 
the Potts model with long-range interactions 

Zvonko Glumac^, 

Faculty of Electrical Engineering, Kneza Trpimira 2B, 31 000 Osijek, Croatia 

Katarina Uzelac 

Institute of Physics, P.O.B. 304, Bijenicka 46, HR-10000 Zagreb, Croatia 
Abstract 

The zeros of the partition function of the ferromagnetic (/-state Potts model with 
long-range interactions in the complex-g plane are studied in the mean-field case, 
while preliminary numerical results are reported for the finite Id chains with power- 
law decaying interactions. In both cases, at any fixed temperature, the zeros lie on 
the arc-shaped contours, which cross the positive real axis at the value for which the 
given temperature is transition temperature. For finite number of spins the positive 
real axis is free of zeros, which approach to it in the thermodynamic limit. The 
convergence exponent of the zero closest to the positive xeal-q axis is found to have 
the same value as the temperature critical exponent 
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1 Introduction 



Several decades ago, it was shown that the study of distribution of the partition- 
function zeros in the complex plane of an appropriate variable, can provide 
relevant informations on phase transitions. In their pioneering work, Lee and 
Yang [1,2] have shown the connection between the existence of phase tran- 
sition and the distribution of partition function zeros of the ferromagnetic 
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Ising model in the complex plane of a symmetry breaking field. Later on, 
[3,4], similar connection was established for the partition function zeros in 
the complex-temperature plane. In both cases, the singularities at the transi- 
tion point are closely related to the way that the zeros approach (either with 
temperature or with the size of the system) to the positive real axis. 

The partition function zeros were investigated for number of models, including 
the Potts model both in complex-field [5,6] and in complex-temperature [7-10] 
plane. In the Potts model a relevant quantity is also the number of states 
often considered as a continuous parameter, which is formally possible on the 
basis of its graphical representation [11] as it will be discussed later in text. 
This gives the motivation to investigate the zeros for this model also in the 
plane of complex q. Such an idea was already formulated [12,13] but mostly 
in context of chromatic polynomials in connection with studies of the ground 
state of the anti-ferromagnetic Potts model [14-20]. Only a few works on 
complex-g plane address the ferromagnetic Potts model at finite temperatures. 
The papers by Chang and Shrock, [21], study the loci of partition function 
zeros in 2d Potts model with different boundary conditions. In their very 
recent paper Kim and Creswick [22] give a systematic study of both Fisher 
and complex-Q' zeros in the 2d Potts model with short-range interactions and 
find the similarity between the scaling property of complex-gr zeros and the 
den Nijs expression for the thermal critical exponent. 

In the present paper we investigate the zeros in the complex-g plane in con- 
nection with the first- and second-order phase transitions in the ferromagnetic 
Potts model. To this purpose we consider here the Potts model with long-range 
power-law decaying interactions, which already in one dimension displays di- 
versity of critical regimes and embraces most of the nontrivial aspects of phase 
transitions, characteristic of the Potts model in general. The hmiting case of 
this model, where all the interactions are equal (and which represents the 
mean-field case of the model) has the advantage to be exactly solvable for 
arbitrary q. 

In the next section we describe briefiy the starting Hamiltonian and its graph- 
ical representation. In subsequent sections we present the results for the mean- 
field (Section 3) and the power- law decaying interactions case (Section 4). 
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2 Hamiltonian 



The ferromagnetic Potts model with interactions of arbitrary range is de- 
scribed by the Hamiltonian 

H^-Y, Ji,j-S{si,Sj), (1) 

where Si denotes the Potts spin at site i of a d-dimensional lattice, which can 
take q values, and Jjj > denotes the ferromagnetic interaction between spins 
at sites i and j. 

An impressive amount of studies of this model in the case of nearest neigh- 
bour interactions (for a review see Wu [23]) reveal its very complex critical 
behaviour which varies with q and comprises transitions of different orders 
and belonging to various universality classes. 

The model with long-range interactions was studied less, and mostly in two 
special cases. The first is the case when all the interactions are equal and 
it represents the mean-field limit of the model. It can be approached by the 
saddle point approximation [24-27] and in the thermodynamic limit the exact 
expressions, in particular those for critical temperature and exponents may be 
derived for all values of q. Second is the case with power-law interactions de- 
caying with distance as The type of phase transition there depends in 
addition on the parameter of range a, which in certain aspects has similar ef- 
fects as the change of dimensionahty [28] . The phase transition is there present 
in one dimension [29-31] and by varying a, series of different regimes are met, 
which makes this Id model a useful paradigm for most of important features 
of the model in general (including also the first-order phase transition). It was 
a subject of number of recent studies [31-39]. 

The graphical representation of the Potts model [11] holds for arbitrary range 
of interactions. The partition function for the Hamiltonian (1) may be ex- 
pressed as 

= E n (2) 

all Q active 
links 

where the summation is taken over all possible graphs Q. Each graph represents 
one possible configuration of connections between the N spins, called active 
links, each of which contributes with a factor Vij = exp[Jij/kBT) — 1. n{Q) 
denotes the number of of disconnected parts of the graph called clusters. The 
product is taken over all the active links of a considered graph. 
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Whatever the interactions are, the number of clusters that a single graph 
contains ranges from one cluster (low temperature limit, where all spins are 
interconnected), up to at most clusters (high temperature limit where no 
active links are present and each spin is a cluster of size one) . Thus, by col- 
lecting all the graphs with the same number of clusters, the partition function 
becomes simply an A'"-th order polynomial in q 

N 

ZN^J2('n{T,Ji,j,N)q^. (3) 

n=l 

This made possible to treat the number of states g as a continuous parameter. 
It also gives the possibility to examine the zeros of the partition function as 
the zeros of the iV-th order polynomial in q. 

The coefficients a„ are functions of temperature, the finite size and interac- 
tions, and are quite complicated in the general case (1), but can be written in 
more compact form for the special cases considered here. 



3 The mean-field case 

The simplest case of long-range interactions is the one where all the interac- 
tions are set to be equal, Jij — J, which reduces the Hamiltonian (1) to the 
form 

N-l N-i 

Hmf = HMF/ksT ^ -— E E ^ (^i' ^i+i)' (4) 

where K = J/ksT, and division by N is necessary to keep the summation 
finite. This model corresponds to the mean-field hmit. It was solved by the 
saddle point approximation [24] . The phase transition is there of the first order 
for q > 2 and of the second order for q < 2. The exact expressions for the 
inverse of the transition temperature Kf is given by 

K, = 2^—^ln{q-l), (5) 
q z 

for g > 2, and by Kc — q ior q < 2. Exact analytical expressions may be derived 
also for the critical exponents matching the values obtained in the mean-field 
regime of the short-range interaction model [27] (except those related to the 
correlation function, as it will be discussed later in text). 
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The graph representation (2) reduces in the mean-field hmit to 



(6) 



axig 



where the active hnks have all the same strength v — exp{K/N) — 1, inde- 
pendently of distance between the spins they connect and b{G) denotes the 
number of active links in a given graph Q. Written as a polynomial in q the 
above equation is equal to 



where the coefficients a„ are real and positive and depend on temperature 
and size N only. (Notice, that, since the number of active bonds may reach 
up to N (N — l)/2, the above partition function could equally be written as a 
polynomial in temperature variable v, which would be of order N {N — l)/2.) 

3.1 Numerical results 

For relatively small sizes the evaluation of coefficients and finding the zeros of 
the polynomial (7) in the complex plane may easily be performed numerically. 
Numerical results for sizes up to N — 35 spins are presented in Figs. 1 (a)-(c) 
for the three different temperatures, K — 0.5, 2, 4.540457. These temperatures 
are the transition temperatures of the mean-field model with q — 0.5, 2, and 
8 respectively. Each figure displays the data obtained for several sizes. 

One may observe that in all the three cases the zeros lie on contours similar 
but different from circle which cross the real axis near the value of q, for which 
the given temperature is the transition temperature. The zeros are absent from 
the positive real axis, but approach to it with increasing N. 

We were thus interested to perform a more precise extrapolation of this cross- 
ing point in the thermodynamic limit and studied the convergence of the zero 
closest to the positive real axis with A^. Further, we have examined the pos- 
sible connection of the convergence of the closest zero with N to the scaling 
properties of the model near criticality. 

To this purpose we have used Burhsch and Stoer (BST) extrapolation pro- 
cedure [40,41]. This algorithm was shown [41] to converge rapidly, needs a 
modest amount of data and, in advantage over the similar algorithm by van 
den Broeck and Schwartz (VBS) [42], it is less sensitive to the rounding errors. 



N 



(7) 



n=l 
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In short, BST algorithm extrapolates a sequence of the form 



A{N) = A^ + ai- iV-" + a2 ■ iV'^- + • • • (8) 
in the N ^ oo limit by a set of transformations recursively given by 



(9) 

, /j(iV+l) _ AN) \ _ 



N J \ 4(^+1) 4(^+1) 



where ^4^^ = and A^^^ — qi^^{N) are the input data and a; is a free param- 
eter. 

If the original series is of the form: 

A{N) ^A^ + ci ■ N-^' + C2 ■ N"^^ + • • • (10) 

it incorporates the leading convergence exponent as a;, while the higher order 
approximants converge as N~'^ , where uj — min{2ui,uj2)- Although it was 
shown, that the extrapolations are insensitive to the choice of a; in a wide 
range of values, it influences its convergence speed. Thus, by maximising the 
convergence of higher approximants one can in the same time obtain both 
the extrapolated value and the leading convergence exponent. Also, when the 
extrapolated value is known in advance, uj yields the leading convergence ex- 
ponent. 

We have examined the convergence of the real and imaginary parts of the zero 
closest to the real axis, denoted as Qi^li^) several temperatures using the 
numerical results ranging from A'^ = 21 to A'^ = 35. 

In the extrapolation of both real and imaginary part qi^l{N), the same pa- 
rameter ou was used and was fixed by the constraint q2'\N — > 00) ^ 0. 

The results of extrapolations are presented on Table 1. By Qo is denoted the 
value of q for which the considered temperature is known to be the tran- 
sition temperature. All the data were obtained by requiring that condition 
g2°'*(A") should be fulfilled with precision of 10~^. 

The presented results indicate that: 

(i) loci of complex-g zeros intersect the positive real-g axis at value qi — Qq- 
The values agree up to a few percents. 
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Table 1 

The BST extrapolation results of the positions of closest zeros. 



K 


0.1 


0.5 


1.0 


2.0 


2.772589 


4.540457 


Qo 


0.1 


0.5 


1.0 


2.0 


3.0 


8.0 


Qi 


0.09998 


0.504 


0.95 


2.05 


3.02 


8.00009 




0.34 


0.33 


0.33 


0.55 


0.91 


1.19 



(ii) Power-law convergence of the closest zero is observed in all the cases 
considered. 

The convergence exponent uj represents the scaling with size of the distance 
from the real axis to the closest zero. Before proceeding with the finite-size 
scaling analysis it is important to notice that here the scaling is not done with 
the linear size, but rather with volume, more precisely, with number of spins, 
since in the infinitely coordinated model the distance as well as the dimension- 
ality lose their meaning. In this respect the considered mean-field limit differs 
from the standard mean-field approach to the short-range interaction model. 
As pointed out by Botet et al. [43,44] the basic scaling quantity in the infinitely 
coordinated system is not the correlation length, but the "correlation number" 
Nci which scales with temperature with the critical exponent v* = vmf ■ dc- 
The exponent umf and dc are, respectively, the mean- field value of the critical 
exponent u and the upper critical dimension (the dimension above which the 
mean-field regime sets on) of the short-range model. 

Thus, the temperature critical exponent for the present model, where only the 
scaling with the number of spins can be applied will be l/u* and not 1/umf- 
In the Potts model vmf — 1/2 both for 5 < 2 and q — 2, while dc — ^ for 
q — 2, but dc = 6 for g < 2 for symmetry reasons. This gives v* — 2 lor q — 2 
and V* — Zior q< 2. 

The results presented in Table 1 clearly show that in the regime of second- 
order phase transition (g < 2) the convergence exponent coincides up to a few 
percents with the temperature critical exponent Xjv* . This result should be 
compared to the one in complex temperature plane, where it was shown [45] 
that the distance of the zero closest to the real axis scales as the tempera- 
ture critical exponent. Similar behaviour in a complex-g plane was recently 
observed in the 2d ferromagnetic Potts model by Kim and Creswick [22] , but 
the precision of their results did not permit them to be conclusive whether 
the obtained convergence exponent was indeed equal to the exponent \jv. 

In cases g = 3 and g = 8 we obtain approximately linear convergence, propor- 
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tional to the system size, which is the form of scahng that is expected for the 
first-order transition. 



3.2 Large N expansion 



In the hmit N ^ 1 some analytical results may easily be derived by extending 
the mean-field solution by Kihara et al. For the large number of spins N , 
the free energy density, without the terms which are constant in the order 
parameter is given by 



T(K,q,S)= (11) 

-S^f^ + ii^^ln[l + ta-l)Sl + fe-l)i^Ml-S). 

z q q q 

The partition function is given by 



1 



Zn{K, q)^ j exp {-N T) dS, (12) 



where the proportionality sign stands for all the omitted constant terms. 

We limit our considerations here to the first-order regime of temperatures, 
where the free energy density at the transition shows two equal minima, and 
the partition function integral can be approximately solved by a saddle point 
method 

Zj,^ l + exp[-NJ^{Smin)]- (13) 



The first term on the right hand side denotes the minimum at 5" = 0, and 
the second term is the contribution of the second minimum at 5" — Sjjii^ > 0, 
which is obtained from the equation for the extremum 

^5.„i±i^. (14) 



The straightforward way to find the complex-g zeros would consist of two 
steps: first, by using (14) one should represent as a function of K and q 
only; second, one should solve the equation Z^r = by analytic continuation 
of q to the complex plane. Since T — J-r + iJ-i then becomes complex, the 
vanishing of Z^, as given by (13), requires 

J^R = 0, NJ^i = ± (2m + l)7r, m = 0, 1, . . . (15) 
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These two equations determine the position of complex conjugate pairs of 
zeros ± ^2"*'*) in the complex-g plane. 

The first of the above mentioned steps cannot be performed in the straight- 
forward way. Instead, we will solve (14) in variable q 

5=^[exp(i^5)-l], (16) 



and represent JF as function of S and K. 



T(K^) ^ (1 - ^) exp(K^) + 1 



The equations (15) are then solved in the complex-S* plane and, by equation 
(16), mapped back into the complex-g plane. 

The solutions gj'g^ should comply with both equations (15). Since the first 
of equations (15) is A'"-independent, it contains all the possible positions of 

zeros for arbitrary N including the limit ^ 00. The zeros of the system 
of particular size N are then determined as intersections of this curve with 
solutions of the second equation in (15). For large these intersections appear 
to be equally distributed along the the curve J-'r = 0, so that it can be drawn 
as the contour containing zeros, (see Fig. 2 as an illustration). The curve 
crosses the positive real-g axis at a point which is transition point for a given 
temperature. 

Since the locus of zeros has a shape similar to a circle, we draw the circular 
contour (long dashed line) in order to stress the difference. The parameters of 
this circle are determined by the facts, that its center hes on the real-q axis 
(since zeros appear in complex conjugate pairs only) and that it intersects the 
positive rcal-gf axis at q solution of (5) and the negative real-q axis ai q — — \q\, 
solution of 

^ = 2M±lin(|g| + l). (18) 



We were not able to give analytical expressions for all the zeros of the model, 
but the location of zeros closest to the positive real S ot q axis is a relatively 
easy task. 
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Since we limit here to the first-order transitions, the temperature will be fixed 
through the parameter Qo > 2 



i^ = 2^^1n(go-l), (19) 



In the thermodynamic limit the second minimum of J-'{S) appears at So — 
{Qo — 2)/{Qo — 1). For large but finite A^, we will expand free energy (17) 
around the minimum, S = So + s. The small s <^ So vanishes as N — > oo, so 
we expect to have, to leading order in 1/N, 

S = Si+lS2, 51 = ^, S2 = ^, (20) 



for real constants a,b,xi,X2- The expansion of (17) and the constraint (15) 
give the equation 



T^iQo- 1) [1 - In (go - 1) (^1 - 5i)] s + e»(s2) = ± i TT, 
( K-Qo + l 2 \ Qo - 1 J. , .... 

The closest zero has the smallest imaginary part, so its index is m = 0. By 
inserting (20) in the above equation, one reads 



while xi is not defined. Mapping back to the complex-g plane by using the 
equation (16) gives 



qf\N)=Qo + 

(Qo- 



+ 



(Qo - 1) 

Qo-2 



[K-Qo) 



Qo-2 



2 (Qo - 1) 



Qo-2 

(Qo - 1)^ 

Qo-2 



Qo) 



NX2 



+ 



^[K-2{Qo 



-11 



jV2x2 j 

(23) 



2ah 
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Finally, by inserting values for a,b,X2 from (22), one obtains the large- A'" 
dependence of the distance of the closest zero from the real axis in the case of 
the first-order transitions. 



^(0)^^) = iK-Qo)^ + , (24) 

The above expressions can be checked numerically. For illustration, the nu- 
merical data of the closest zeros, for the temperature critical for g = 8 are 
presented in Table 2 for different sizes N. 

Table 2 

The location of the closest zeros for the temperature given by the relation (19) for 
Qo = 8. N denotes the number of the system considered. 



N 










5000 


7.99999543669504 


-8.679091351481774 


■ 10- 


-3 


6000 


7.99999683103614 


-7.232577480179549 ■ 


■ 10" 


-3 


7000 


7.99999767178073 


-6.199352825621011 


■ 10" 


-3 


8000 


7.99999821745666 


-5.424434119811201 


■ 10" 


-3 


9000 


7.99999859157044 


-4.821719459791061 


■ 10" 


-3 


10000 


7.99999885917194 


-4.339547669719317 


■ 10- 


-3 


11000 


7.99999905716679 


-3.945043440975050 


■ 10- 


-3 


12000 


7.99999920775815 


-3.616289894012384 


■ 10- 


-3 


13000 


7.99999932495365 


-3.338113800831173 


• 10" 


-3 


14000 


7.99999941794470 


-3.099677139486969 


• 10- 


-3 


15000 


7.99999949296515 


-2.893032025972416 


■ 10" 


-3 



The data q'iJ(A^) from Table 2 and analogous series of data for different tem- 
peratures corresponding to q = 3 and q = 2.1 were extrapolated to the ther- 
modynamic limit by using the BST algorithm. 

The extrapolation parameter oo was fixed by the constraint q2'\N) 0. After 
we have made the extrapolations of qi^li^) ^-iid found the values of ^1(00) 
and LU, we have performed alternative calculation of convergence exponents of 



11 



qi^^(N). By assuming the convergence of the gi^2(-^) to be of the form 

q?\N) - gi(oo) = Ci • 7V-"^ + • • • , qi°\N) = C2 • iV""^ + • • • ,(25) 



the convergence exponents ai,2(-^) for real and imaginary part of the closest 
zeros, can be expressed as 



a2{N) = 



In { [q?\N) - q,{oo)] / [qf'\N + 1) - gi(oo)] } 



\n[N/{N + 1)] 



In 



Q2 



\n[N/(N + l)] 



(26) 



The above data for q;i,2(^) are extrapolated again by the BST algorithm. The 
complete set of extrapolated results is summarised in the Table 3. 

Table 3 

The results of the BST extrapolations performed on the closest zeros, ^1,2 (00), 
followed by parameter u and convergence exponents ai^2, for set of temperatures 
determined by Qo and (19). The predicted convergence exponents xi^2 as given by 
(24) are in the last two columns. 



Qo 


gi(oo) 


92(00) 


LO 


ai 


Q!2 


Xl 


X2 


8 


8.0000000000001 


3 • 10-13 


1.0000 


1.9999988 


1.000000002 


2 


1 


3 


3.000000000005 


4 • 10-12 


1.00000 


2.00003 


0.99999997 


2 


1 


2.1 


2.1002 


9 • 10^6 


0.98 


1.85 


0.95 


2 


1 



We observe a full agreement between the analytically given values for q and 
convergence exponents on one side, and the numerically calculated and extrap- 
olated values for q and exponents on the other side. Although the scaling in 
the case of the first-order phase transition is trivial and linearly proportional 
to system size, the results can be used to check the convergence of numerical 
results. For g = 3 and q = 8 the convergence exponent has improved compared 
to the earlier numerical results given in Table 1, confirming that larger error 
in this case comes from the small size of considered systems and that scahng 
is indeed linear with size as expected for the first-order phase transition. The 
loss of precision in extrapolations persists only for the case corresponding to 
q = 2.1 and should be attributed to the crossover effects due to vicinity of the 
border between the first- and second-order transitions. 
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For the partition function zeros lying on the curve in the appropriate complex 
plane, the normalised density of zeros g, is defined as [2] 



where dn is the number of zeros inside the arc of the curve of the length 
dl. In the case of the field- or temperature-driven second-order transitions, 
the density of zeros vanishes at the transition with an exponent connected 
to critical exponents [46]. On the contrary, at the first-order transition point, 
the density of zeros remains constant. To the end of this paragraph, we will 
examine the behaviour of density of complex-g zeros in the vicinity of first- 
order transition point. To this purpose, we notice that the large- expansion 
of qi,2{N) for the closest zeros (those with m — 0) presented earlier in text, 
may be directly generalised to other values of m as long as the condition 
(2m + 1) < is fulfilled. This leads to 



gi \N)-q^{cG)^C^ • ( ^ j +•••, 

2m + 

N ' ^ 



(28) 



with the constants gi(oo),Ci and C2 given by the relations (24), where 1/N 
is changed by (2m -|- 1)/-/V. 



To calculate numerically density of zeros (27), we choose two neighbouring 
zeros, with dl denoting geometrical distance between them. 



dl 



(m+1) (m) 
?1 - ?! 



+ 



(m+1) (m) 
H2 ~ H2 



near positive real-g axis becomes 



In that case density of zeros 



, , 1 [^2 /2m^ 
9{m) -J^\CI 



+ cl 



N 
N 



1 + 



1 + 



2m + 1 
2 



- 1 



2m -M 



.2^-1/2 

- 1 



(29) 



The edge (m = 0) density of zeros qq, tends to constant value in the thermo- 
dynamic limit. 



^0 = 



2C, 



Cl 1 



2a' 



(30) 
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For example, for Qo = 8 and K{Qo = 8) = 4.540457014462 this constant is 



2 Co 



go-2 



{K-Qo)b 



0.011521939 



(31) 



For illustration, we check this result numerically by calculating for the values 
of Qo and K from above, and N — 10 000, the twenty zeros closest to the 
positive real-g axis and corresponding densities. They are presented in Table 
4. 

The limiting value of the density at the transition point, is extracted by the 
BST extrapolation procedure (8) simply by changing variable 1/N into (2m + 
1)/N. 

Qm^go + g [ ) +•■■■ (32) 



The extrapolations give go = 1.152194 ■ 10 ^ for a wide range of 0.3 < a; < 3.0, 
which reproduces the analytically obtained value (31) in six significant digits. 

We can conclude that the density of the partition function zeros in the complex- 
q plane behave in similar way as the density of zeros in the complex-temperature 
[3] or complex-field [46] planes: in thermodynamic limit it has constant value 
at the first-order transition point. 



4 The power-law decaying interactions 



The second case that we have considered is the ferromagnetic Potts model in 
one dimension with power-law decaying interactions. When taking the periodic 
boundary conditions, it is described by the reduced Hamiltonian 



H 



PL 



N-l N-i 



+ 



(33) 



Although in one dimension, the model is nontrivial and has a phase transition 
at nonzero temperature for all q when < cr < 1 [29-31]. This transition is of 
the mean- field type for low enough values of cr, cr < crc{q), and is a non-trivial 
second-order phase transition for cr > crc(Q'). In the classical regime cr < a^q), 
the transition is of the first order when q > 2 [34,35,37]. The exact position of 
the line separating the two regimes, crc(g), is a difficult and still open question 
[35,37,47]. 
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Tabic 4 

The twenty zeros closest to the positive real-gr axis and the corresponding densities 
are shown at temperature determined from the relation (19) and Qq = 8 and for 
the AT = 10000 spins. 



(m) 

m ql 

7.99999885917192 

1 7.99998973259138 

2 7.99997147966529 

3 7.99994410086352 

4 7.99990759689083 

5 7.99986196868649 

6 7.99980721742419 

7 7.99974334451167 

8 7.99967035159032 

9 7.99958824053479 

10 7.99949701345249 

11 7.99939667268302 

12 7.99928722079755 

13 7.99916866059815 

14 7.99904099511710 

15 7.99890422761603 

16 7.99875836158505 

17 7.99860340074195 

18 7.99843934903115 

19 7.99826621062264 



-0.004339547669719 
-0.013018627057442 
-0.021697658591626 
-0.030376610375500 
-0.039055450520507 
-0.047734147149695 
-0.056412668400979 
-0.065090982430549 
-0.073769057415945 
-0.082446861559540 
-0.091124363091702 
-0.099801530274076 
-0.108478331402754 
-0.117154734811648 
-0.125830708875537 
-0.134506222013369 
-0.143181242691353 
-0.151855739426127 
-0.160529680787899 
-0.169203035403569 



g{m) 

0.0115219529910602 
0.0115219974088187 
0.0115220714335465 
0.0115221750581617 
0.0115223082726076 
0.0115224710640261 
0.0115226634165474 
0.0115228853117394 
0.0115231367280630 
0.0115234176412815 
0.0115237280243068 
0.0115240678473335 
0.0115244370775645 
0.0115248356796385 
0.0115252636152771 
0.0115257208435568 
0.0115262073207999 
0.0115267230005962 
0.0115272678338162 
0.0115278417687404 



In the graph representation, the difference from the mean-field case is that 
the contribution of active hnks depends on distance between the spins they 
connect. The contribution of active bonds may then have N/2 (for N even) 
or (A^ — l)/2 (for odd) distinct values given by Vj = exp[fC/j^+°" + K/ {N — 
JY^'^] — 1. If we denote by bj{Q) the number of active links of the length j or 
— j in a graph Q, the partition function becomes 

Z^^Y: ^1^'^ ^2^'^ ■ ■ ■ <f^^^ (34) 

all0 
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or, in form of polynomial in q, 

N 

ZN^J2<'n{vj,N)q-. (35) 

n=l 

The real and positive coefficients a„ depend on temperature and N, but also 
on the parameter of range a. This makes the calculations more demanding for 
the computer memory. 

We have performed exact numerical calculations on chains with < 9 spins 
for two sets of parameters, {a = 0.3, K = 0.576) and [a = 0.8 K = 0.8230). 
As in the previous case, the chosen temperatures are the critical temperatures 
known for certain values of q. They correspond to g = 5 and q = 2 respectively. 
Remark, however, that in the present case the critical temperatures are not 
known exactly. The values used here were obtained by the finite-range scaling 
approach [31]. The ffi^st and second set correspond respectively to the to the 
first- and second-order phase transition. 

The numerical results for the loci of complex-g' zeros for sizes N — 6, - ■ ■ ,9 are 
presented in Figures 3 and 4. In both examples are obtained the zeros lying on 
the arc-shaped curves, while the positive real axis is free of zeros. The zeros 
closest to it approach to it with increasing N. 

As in the mean-field case, we have examined the convergence of the closest 
zeros to the thermodynamic limit by the BST algorithm, assuming the power- 
law form. The analysis included the data obtained for chains with = 5 to 
9. The parameter oo was determined, as earlier, by requiring that g2 ^ in 
the thermodynamic limit. The results are presented in Table 5. 

In the case of second order phase transition the expected crossing point with 
the real axis is obtained with precision better than 1%. The leading conver- 
gence exponent given by uj, is equal to 0.41, which is close to the value for 
the temperature critical exponent l/v = 0.48 obtained earlier by finite-range 
scaling for q = 2 and a = 0.8 [31]. 

In the case of the ffist-order transition the fit is more difficult. The discrepancy 
between the obtained crossing point and the value expected for the considered 
temperature is larger, about 20%, while the convergence exponent, which in 
context of finite-size scaling should be linear for a first-order phase transition, 
differs by about 11%. Larger error in determination of the value of q, may be 
attributed to both lower precision of the data for critical temperature obtained 
by finite-range scaling [31,35] when the values of a are low and to the proximity 
of the second-order transition regime in the (g, cr) plane, which might induce 
crossover effects for small sizes. 
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Table 5 

The loci of closest zeros and its BST extrapolation at {K = 0.576, a = 0.3) and 
{K = 0.823, a = 0.8) 





K = 0.576 


a = 0.3 




N 


HI v'y 








-f) 44-741 '^89QnQ'^9Q8 






Q 


0734103695976816 


2 88316232503579 




7 


0.486108948338677 


2.88942370075814 




8 


0.818954244961081 


2.85269647435491 




9 


1.09269817849075 


2.79489503015230 




oo 


4.2 


3 • 10"^ 


0.89 




K = 0.823 


a = 0.8 




N 

1 V 








5 


-0 418572384698087 


3 69577294340190 




6 


0.179124021423254 


3.55212427658784 




7 


0.581821465484680 


3.37723067427777 




8 


0.866098195024489 


3.20669717997192 




9 


1.074705295884560 


3.05042285906340 




oo 


2.01 


3 ■ 10-=^ 


0.41 



Since in these preliminary calculations we dispose with small number of data, 
we have not calculated the density of zeros in the case of power-law model. 



5 Conclusion 

We have studied the complex-g zeros of the partition function of the ferro- 
magnetic Potts model with long-range interactions in context of the first- 
and second-order phase transitions in this model. Two cases are presented: 
the mean-field limit and a few examples of the case with power-law decaying 
interactions. 

We conclude that the positions of the zeros may be analysed in a similar way 
as the zeros in the plane of complex field or temperature and that they can give 
useful informations on the phase transition in ferromagnetic Potts models. 

Numerical results for relatively small systems show the similar picture as in 
the case of zeros in the complex field or temperature plane. In all the cases 



17 



considered the zeros lie on smooth curves that intersect the real positive axis 
at value for which the given temperature is the transition temperature. For 
finite systems the positive real axis is free of zeros, but they approach to it 
with increasing size. 

The distance of the closest zero from the real axis is found to scale with size 
as the temperature critical exponent 1/v. In the mean-field case, this result 
could be obtained with a good precision for series of different temperatures 
corresponding to different transition regimes. In the case of power-law inter- 
actions, the preliminary results obtained with less precision suggest the same 
behaviour. This result is similar to behaviour in the complex temperature 
plane, and it supports recent findings by Kim and Creswick [22] in the short- 
range Potts model in 2d. 

In the end, let us mention some technical advantages of the analysis by 
complex-g zeros compared to the complex-temperature zeros in the model 
with long-range interactions. The partition function of the mean-field model 
consisting of spins can be written as a iV-th order polynomial in g, or 
roughly iV^-th order polynomial in temperature variable. Consequently, this 
allows better numerical precision of calculations in the complex-g plane then 
those in the complex-temperature plane. For a model with power-law decay- 
ing interactions the possibilities in complex temperature plane are even more 
restricted. Model with spins, has N/2 (for N even) different parameters 
involving temperature, so that the partition function cannot be defined as a 
polynomial in some temperature variable at all. 

Further investigations of the complex-q' zeros should be done in future in at 
least two directions. We plan to perform a more extended numerical calcula- 
tions of the complex-g zeros in the case of power-law decaying interactions. 
This should allow us to obtain the convergence exponent for the closest zero 
with more precision. In addition, it would be interesting to examine a still 
open question of the borderline qdcr) dividing the first- from the second-order 
transition regimes within the formalism of complex-g' zeros. Also, on the basis 
of the results obtained in this paper, a possible connection between the zeros 
in complex-g and complex-temperature plane deserves further investigation. 
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Figure caption 



Fig. 1: The loci of complex-gf zeros of the partition function of the mean-field 

Potts model calculated at inverse temperatures (a) K = Kc{q = 0.5) = 0.5, 
(b) K = KrXq = 2) = 2 and (c) K = Kt{q = 8.0) = 4.540457 for = 25 
(open diamonds), = 30 (filled triangles) and = 35 (filled circles) spins. 

Fig. 2: The locus of complex-g zeros of the partition function of the mean-field 

Potts model calculated at inverse temperature K = Kmf{q = 8) = 4.540457. 
The full line denotes loci of zeros in the limit N ^ oo, obtained by the 
approach by the saddle point approach as a solution of the equation jF^j = 0. 
Long dashed line is a circle (described in text) drawn as a guide to the eye. 

Fig. 3: The complcx-g zeros of the partition ftmction of the Id Potts model 
with power-law interactions at the FRS estimate of first-order transition tem- 
perature K = Kt{q = 5, (T = 0.3) = 0.576, for N = 6 (open diamonds), N = 7 
(filled triangles), N = 8 (open squares) and N = 9 (filled circles) spins. 

Fig. 4: The complex-g zeros of the partition function of the Id Potts model 
with power-law interactions at the FRS estimate of second-order transition 
temperature K = Kc{q = 2,0" = 0.8) = 0.8230, for = 6 (open diamonds), 
N = 7 (filled triangles), N = 8 (open squares) and A^ = 9 (filled circles) spins. 
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Fig. 1 (a) 
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Fig. 1 (b) 
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Fig. 1 (c) 
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Fig. 2 
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Fig. 3 
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Fig. 4 
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